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Abstract 

The current-voltage characteristics of Superconducting Quantum Interference 
Devices (SQUIDs) are known to modulate as a function of applied magnetic 
field with a period of one flux quantum <I>o = h/2e. Here we report on the 
fabrication and properties of SQUIDs modulating with a period of 1/2 x <I>o- 
The characteristics of these bicrystal SQUIDs are consistent with either a 
strong sin(2(/j) component of the current-phase relation of the Josephson cur- 
rent, or with an interaction between the Cooper-pairs, causing an admixture 
of quartets to the condensate. 
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The interference between two superconducting condensates results in a wide range of 
intriguing, yet easily accessible phenomena. As shown by B.D. Josephson [1], interference 
between the superconducting order parameters controls the current- voltage {I{V)) char- 
acteristic of weak links. In first approximation, the density of the zero-voltage current 
traversing such a link is given by: 

J = JoSm{Aip), (1) 

where Jc is the critical current density of the junction, and Acp describes the difference of 
the phases of the superconducting order parameter on both sides of the junction. In case a 
voltage is generated by the junction, it is proportional to dfAip : 

dAif = e*vin, (2) 

where e* = 2e is the charge of the charge carriers in the condensate. Embedding one or two 
Josephson junctions into a superconducting loop that encloses an area A, superconducting 
quantum interference devices (SQUIDs) are obtained [2,3]. Exploiting the phase shifts 

= {e*/h)A, (3) 

induced by the magnetic flux density hqH penetrating the SQUID loop, these devices are 

— * 

magnetometers with outstanding sensitivity. Here, A is the vector potential associated with 

— * 

H. According to Eqs. 1 and 3, the circulating currents in the SQUIDs, and thereby the 
SQUIDs' characteristics, periodically modulate with H, the period A$ being given by the 
magnetic flux quantum $o 

A$ = $0 = h/e* =/i/2e. (4) 

SQUIDs have been fabricated in huge numbers. In all cases in which the oscillation 
period was measured, it was found to equal h/2e. In fact, rings containing grain boundary 
Josephson junctions have been used with great success to measure the Cooper-pair charge 
in YBa2Cu307_5 [4]. Therefore one might believe that the oscillation period of SQUIDs 
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always has to be $o- This behef would be erroneous, because neither is A$ — h/2e required 
for SQUIDs that use junctions for which Eq. 1 is not fulfilled (see, e.g., [5]), nor do Eq. 1 
and 3 predict $o-periods for SQUIDs built from hypothetical superconductors with carriers 
characterized by charges e* ^ 2e. 

Intriguingly, experiments measuring the impedance of rf-SQUIDs [6] and transport mea- 
surements of dc-SQUIDs [7] indicate the contribution of a sin(2(/7) component to the Joseph- 
son current for (001)/(110) grain boundaries in YBa2Cu307_5, such that for these junctions 
the relation between the Josephson current density and the phase difference is given by 
J = Jcisin(A(/7) — Jc2 sin(2A(/7). Interestingly, for faceted (100)/(110) boundaries it was 
noted that the spatially averaged current density may consist of many harmonics, even 
if the local tunnel current density only has a first harmonic component Jc2 [8]. In case 
Jc2 dominates J, SQUIDs built from such junctions are expected to reveal a periodicity 
A$ = 1/2 X $0- 

The second possibility to generate non-/i/2e periods, superconducting carriers with 
charges e* ^ 2e, seems to be implausible. Nevertheless, superconductors with such car- 
riers can in principle exist. We were led to consider this idea while measuring the magnetic 
flux generated by (100)/(110) boundaries in bicrystaUine YBa2Cu307_5 films [9]. As shown 
by Fig. 1, in one of these samples the flux is generated in quantities close to $ ~ /i/4e and 
$ ~ /i/6e. Although these data could not be confirmed with other samples, we analyzed the 
possibility that in superconductors such as the cuprates an interaction between the Cooper- 
pairs binds a fraction of the pairs into particles with charge e* = nx2e, n = 2,3,4... . 

These considerations suggested the fabrication of SQUIDs from (100)/(110) boundaries 
to analyze their oscillation period. The need for such measurements is underlined by mea- 
surements of K. Char etai., who found Shapiro steps [10,11] at (100)/(110) boundaries to 
occur not only at the standard voltages Vn — nh/2ef, but also at voltages Ki = nh/Aef [12]. 
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Here, / is the microwave frequency used to generate the resonance steps. Furthermore, 
studying the critical current /c of such a boundary in a Bi2Sr2CaiCu208+a; film, D. van 
Harlingen and his group noted that the width of the Ic{H) zero-field maximum had only 
half the expected value [13]. 

While it is possible to attribute these experimental results to effects well known in su- 
perconductivity and to peculiar artifacts resulting from inhomogeneities, we saw the need 
to analyze the oscillation periods of dc-SQUIDs built from asymmetric 45° boundaries to 
clarify whether they equal h/2e. We note that numerous SQUIDs containing such junctions 
have already been fabricated with the so-called biepitaxial process. To our knowledge a 
special oscillation period was not reported for any (see, e.g. [14]). This behavior may be the 
result of Jc inhomogeneities masking the true oscillation period. SQUIDs with asymmetric 
45° junctions have also been built using the more reproducible bicrystal technology. The 
dynamics of these SQUIDs demonstrated the presence of strong second harmonics in the 
current-phase relations of the junctions [7]. 

We fabricated and analyzed four bicrystalline dc-SQUIDs containing (100)/(110) bound- 
aries, together with reference samples for calibration. The experiments were performed with 
YBa2Cu307_5 films (3 samples) and with a Yo.6Cao.4Ba2Cu307_5 film (1 sample) grown 
by pulsed laser deposition to a thickness of ~120nm. The commercial SrTiOs-substrates 
contained (100)/(110)-asymmetric grain boundaries, specified to within 1°. The films were 
patterned into dc SQUID structures by photolithography and ion-beam etching before Au 
contacts were sputter deposited. As reference samples, nominally identical SQUIDs were 
grown on symmetric 45° bicrystals. The junctions have widths of 9 /xm, the SQUID holes 
are symmetric triangles with base lengths of 8 //m and heights of 9.5 /xm. A sketch and a 
micrograph of the devices are given in Fig. 2. All measurements were conducted in a magnet- 
ically shielded room. The SQUIDs had non-hysteretic current-voltage /(V^)-charactcristics, 
their critical currents were measured by tracing the /(\^)-curves, using a voltage criterion of 
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~ 2 jjM. Because the Yo.6Cao.4Ba2Cu307_5 sample behaved identically to the YBa2Cu307_5 
SQUIDs, we do not differentiate the two in the following. 

The grain boundaries in the bicrystalline YBa2Cu307_5 films are faceted, with facet 
lengths of the order of 10 — 100 nm. Together with the (\x'^-y2 symmetry of the YBa2Cu307_5, 
the facets cause anomalous Ic{H) dependencies [15]. These arise because a large fraction of 
the facets act as 7r-Josephson junctions with negative Ic- 

Fig. 3 displays the Ic{H) dependence of a SQUID with an asymmetric 45° boundary mea- 
sured at 4.2 K (a), together with the corresponding characteristic of a control SQUID with 
a symmetric 45° boundary (c). Both curves show the SQUID oscillations superimposed on 
the Ic{H) curves of the Josephson junctions. These interferences occur because the width of 
the Josephson junctions is comparable to a side length of the SQUID hole. It is noted that 
for the SQUID with the asymmetric boundary, 7c is small at small H, which is due to the 
suppression of Ic by the vr-facets. /c(4.2K) ai H = typically equals 10 — 100 A for the 
YBa2Cu307_5 samples, which corresponds to Jc(4.2K) = 10^ — lO^A/cm^, typical values 
for 45° boundaries [15]. In these respects, the SQUIDs behave as expected. But how large 
are the oscillation periods? 

As shown by Fig. 3, the SQUID oscillations of the control sample have a period of 
~ 2.7 iiT. To evaluate the magnetic flux in the loop that corresponds to this applied fleld, we 
cahbrated the flux focussing of the samples. To do so, 24° and 36° bicrystaUine YBa2Cu307_5 
fllms were structured into the sample pattern. Their oscillation periods are 2.6 — 2.9 //T [16], 
very similar to that of the control sample. Because bicrystal SQUIDs with almost identical 
geometries have been reported to oscillate with h/2e [17], it is evident that the oscillation 
period of the control SQUIDs is h/2e, corresponding to a flux focussing factor of 20.2. 

In contrast to the standard SQUIDs, the oscillation periods of the SQUIDs with the 
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(100)/(110) boundaries amount to 1.44 //T (Fig. 3a,b). Because the flux focussing factor 
and the loop areas are nominally identical to those of the control SQUID, this periodicity 
equals $o/2 = /i/4e. This period we observed in three out of the four SQUIDs, but never 
in a control sample. In the SQUID shown in Fig. 3a, the period even changes from /i/4e 
to /i/2e, when H exceeds 5 /xT (Fig. 3). Providing an internal cahbration, this curve gives 
additional evidence that the oscillation period equals $o/2 at small H. Both oscillation 
periods are clearly revealed by the Fourier-transform of the /c(-f^)-dependence (Fig. 4). The 
spectra even suggest minute signals at h/6e and h/8e, although the higher harmonics die 
off quickly. The temperature dependencies of both oscillation periods differ. While at 77 K 
the $o/2 have almost completely disappeared, the $o oscillations are still pronounced (see 
Fig. 3d). 

Clearly the dc-properties of the SQUIDs are highly unusual and do not agree with 
simple expectations based on the dc-Josephson relation (Eq. 1). To analyze whether the 
rf-properties of the SQUIDs also differ from the standard behavior as predicted by the ac- 
Josephson-relation (Eq. 2), we measured the I{V) curves of the samples under microwave 
radiation at 11.88 GHz. Since the radiation caused the samples to become noisy at small 
H, these experiments could only be performed in the large H regime, for which the sam- 
ples showed the h/2e period. The I{V) characteristics display conventional, integer Shapiro 
steps Vn — nh/2ef in case the samples are biased with magnetic fields for which the SQUID 
oscillation have /p-maxima. Shapiro steps at half the standard frequency Vn = nh/2ef are 
seen if the junctions are operated in magnetic fields that result in /p-minima. 

The SQUID properties strongly suggest that they are caused by a dominating Jc2- 
contribution to the Josephson current, in agreement with Ref. [6,7]. For (100)/(110) bound- 
aries a large Jc2/ Jci ratio is expected, because Jd is suppressed by the tunneling geometry, 
which couples the nodal with the anti-nodal directions of the order parameters. Further, 
because V<^ is insignificant at small H (Eq. 3), the 7r-facets have negative current densities 
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in this field range: Ji = Jd sin( Ac/? + tt) = — Jd sin(A(/7) and therefore suppress Id- How- 
ever, these facets do not suppress /c2, because at the 7r-facets J2 = Jc2sin(2 x (Acp + tt)). 
Therefore, the Ic2/ hi ratio is enhanced at small i7, precisely where the $o/2-period is seen. 
The suppression of Id by the asymmetric 45° boundaries also explains why the /i/4e-period 
has only been seen for SQUIDs fabricated from such junctions. Because the facet structure 
of the grain boundaries plays a key role in this mechanism, it is understandable that the 
fourth sample we have studied showed A$ = /i/2e, as did the SQUIDs fabricated by Lind- 
strom and coworkers. Due to non-identical growth parameters, the devices vary in their 
facet structures and therefore have different Jc2{x) / Jci{x) ratios. 

Of course, if the data were interpreted in conventional terms, such that the SQUID os- 
cillations are taken as a measure of e*, the /i/4e-oscillations would propose that carriers 
with charge 4e are present at the grain boundary, and that their condensate couples the 
phases of the two Josephson junctions around the SQUID-hole. The possibihty that this 
model was correct seems to be remote, because the BCS ground state is composed solely of 
Cooper-pairs [18]. However, the derivation of this ground state is based on the assumption 
that all interactions between the charge carriers are two-particle interactions, described by 
the term Vj^g. This mean-field approximation is clearly appropriate for conventional su- 
perconductors, in which the coherence length ^ is much larger than the spacing d between 
Cooper-pairs. However, because in the cuprates ^ and d are of the same order, and because 
correlation effects are predominant, it is not clear why this approximation should be appli- 
cable to the high-Tc compounds [19,20]. 

The possible attractive or repulsive interaction between an electron or a hole of one 
Cooper-pair with an electron or a hole of a second pair, arising for example from magnetic 
interactions. Coulomb interactions or lattice distortions, results in a pair-pair interaction 
energy Vcp,cp- Depending on the sign of Vcp,cP) the two Cooper-pairs are correlated or 
anticorrelated. For a strongly attractive Vcp,cp interaction they pair. For weak Vcp,cp, 
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these pairs dynamically rearrange themselves in a fluctuating manner. The correlations 
between the pairs result in a non-BCS ground state, expressed by the Ansatz 

* = Q;2^'2e + Q;4*4e + Q;6*6e + -, (5) 

where Q;2^2e is the Cooper-pair term and the higher order terms describe the Cooper-pair 
multiples, such as quartets or sextets. 

We note that an analysis of the temperature dependencies of the penetration depths 
and of the fluctuation contributions to the conductivity and to the magnetization led to 
the proposal that the superconducting phase transition of the cuprates is due to interacting 
pairs [20]. Intriguingly, a related situation exists for the binding of nucleons in light nuclei. 
While pairs of nucleons condense into a superconducting condensate, they may bind by their 
higher order correlations into multiples, in particular into alpha particle-like quartets [21,22]. 
The onset of "quartetting" in nuclear matter has been considered [23] and found to predom- 
inate in the low density regime. 

Because the Josephson current crossing weak links also couples the higher order multi- 
ples, SQUIDs built from superconductors in which |Q!4p/|Q;2p is non-negligible are expected 
to show an /i/4e contribution in their oscillations. For this to occur, the Q;4^'4e term has 
to be sufficiently large to provide phase couphng around the SQUID loop. Although this 
exotic mechanism is consistent with the measured /i/4e-behavior of the SQUIDs, it has to 
provide answers for the questions why /i/4e periods have never been observed for SQUIDs 
with more conventional junction geometries, and why the Abrikosov vortices carry a flux of 
h/2e. To account for these effects, one would need to conclude that |Q;4p/|a2p is enhanced 
at (100)7(110) interfaces. 

One might expect to find non-/i/2e oscillation periods or non-/i/2e vortices also in other 
systems and under different circumstances, for example half flux quanta in the bulk, in par- 
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ticular if further mechanisms exist that cause the unusual SQUID properties reported. 

In summary, we have fabricated three high-Tc SQUIDs which display Ic{H) oscillations 
with a period of 1/2 x h/2e. Two mechanisms are found to be consistent with these char- 
acteristics. The more conventional is based on Josephson currents with a. J — JcSm{2Aip) 
current-phase relation. The other mechanism is based on interacting pairs that form quar- 
tets with charge 4e. Further experiments are required to clarify whether one of these effects 
underlies the SQUID behavior, and whether these mechanisms are possibly realized in other 
systems. 

The authors gratefully acknowledge very helpful discussions with Y. Barash, B. Cheska, 
U. Eckern, A. Faessler, G. Ingold, V. Korting, C. Laschinger, R. G. Mints, H. Rogalla, 
D. VoUhardt, and C.C. Tsuei. This work was supported by the BMBF (EKM-project 
13N6918), by the Deutsche Forschungsgemeinschaft (SFB 484), and by the ESF (PiShift). 
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FIG. 1. Magnetic flux <I> generated at 4.2 K by a (110)/(100) grain boundary in a 
YBa2Cu307_5-film, measured by scanning SQUID microscopy. 
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FIG. 2. Micrograph (left) and sketch (right) of a YBa2Cu307_5 SQUID used in the study. The 
grain boundary is marked by the arrows labelled 'GB'. The blow-up on the right sketches the facet 
structure and the corresponding Josephson current crossing the boundary. 
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FIG. 3. (a) Magnetic field dependence of the critical current of a YBa2Cu307_5 SQUID with 
an asymmetric 45° grain boundary measured at 4.2 K; (b) same data as (a) with the oscillation 
amplitude numerically normalized; (c) field dependence of the critical current of a control SQUID 
with a symmetric 45° boundary; (d) same as (a), but measured at 77 K. 
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FIG. 4. Fourier transform of the SQUID-oscillations shown in Fig. 3(a). 
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